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Abstract 

Based on the properties of the critical and the actual effective masses of sigma and omega 
mesons, criteria to estimate the values of the isoscalar nonlinear terms of the standard relativistic 
mean field model that reproduce stable equations of state in respect to particle hole excitation 
at high densities are derived. The relation between nuclear matter stability and the symmetric 
nuclear matter properties are shown. The criteria are used to analyze in a more systematic way the 
high-density longitudinal and transverse instabilities of some parameter sets of relativistic mean 
field models. The critical role of the vector and vector-scalar nonlinear terms is also discussed 
quantitatively. 

PACS numbers: 21.30.Fe, 21.65.+f, 21.60.-n 
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I. INTRODUCTION 



Relativistic mean field (RMF) models have been quite successful in providing a micro- 
scopic description of many ground-state properties from medium to heavy nuclei. The 
parameters of each model are obtained by adjusting its parameters to a few ground-state 
properties of a set of magic and semi-magic nuclei. The differences between one parameter 
set and the others are not only due to differences in choosing observables, strategies, and 
constraints for the parametrization, but also to the difference in the nonlinear ansatz used. 
Even in recent years, there have still been some efforts to generate new RMF parameter 
sets [l, 2, 3, 4|. The aim is to improve and to extend the applicability of the model for the 
description not only of ground states but also of excited states of finite nuclei and nuclear 
matter. Indeed, even the boundaries are not too clear, the range of applicability of RMF 
models is limited. Therefore, applying the models in regions which are outside their fitting 
window by extrapolation needs special care. 

Recently, considerable attention has been paid in constraining the equation of state (EOS) 
of matter beyond the nuclear matter saturation density (pq) using flow data from heavy-ion 
collision experiments and astrophysical observatories 5|, [g, l7| . In order to be consistent with 
data, the EOS should be relatively soft at moderate densities ( < 5po) and stiff at high 
densities ( > 5po)- Many models, including some parameter sets of RMF models have been 
checked. The author of Ref. [3], also used some model independent calculations of pure 
neutron matter as an additional theoretical constraint to test the validaty of two parameter 
sets of RMF models (NL3 and FSUGold) at high densities. It is found that FSUGold is 
consistent with all constraints except for a high density EOS that appears mildly softer 



than required by astronomical observations 



7|. It seems that some parameter sets of RMF 



models are unable to fulfil the above requirements 5|, |6|, |7|]. At this point, we have sufficient 
reason to re-investigate the parameterizations of RMF models with high density applications 
in mind. 

We studied the high density instabilities of some representative parameter sets of RMF 
models with respect to density fluctuation by observing their longitudinal and transverse 
particle- hole excitation modes Sj. We have found that certain parameter sets are unstable at 
high densities but for some parameter sets with additional vector and vector-scalar nonlinear 
terms, the onset of instabilities can be pushed into a region with quite large matter density 
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(Pb) and perturbed momenta (g). It means that the longitudinal and transverse parts of 
the particle-hole excitation modes of relativistic mean field models depend sensitively on the 
isoscalar nonlinear terms used. Therefore by adjusting these nonlinear terms, the instabilities 
at high densities of RMF models can be avoided. 

In this work, we extend the results of our previous letter Sj by deriving criteria to 
estimate the values of isoscalar nonlinear terms of the standard RMF model (RMF model 
with minimal nonlinear terms) that produce a stable EOS at high densities. We then 
employ these criteria to systematically study the high density longitudinal and transverse 
instabilities of some representative RMF models. In this way, the actual role of vector and 
scalar-vector nonlinear terms can be revealed in a more quantitative manner. These analyses 
can provide a practical hint that should be useful for RMF parameterization. The extension 
to the isovector sector is also important. The question whether the parameter set of the 
standard RMF model with stable nuclear matter EOS at high densities is consistent with 
the above constraints and simultaneously has acceptable predictions for finite nuclei, would 
also be very interesting to investigate, but these points need additional considerations that 
are beyond the scope of this work and will left to a future paper. 

This paper is organized as follow: in Sec. [TTl it is discussed how to estimate the isoscalar 
nonlinear parameter sets of the standard RMF model; in Sec. IIIII we discuss the stability of 
some RMF models; and finally in Sec. llVl the conclusions are drawn. 

II. CONSTRAINING NONLINEAR PARAMETERS OF STANDARD RELA- 
TIVISTIC MEAN FIELD MODELS 

The instability of nuclear matter at high density is generated by particle-hole excitation 
modes. These modes are indicated by the existence of poles of the meson propagators at zero 
energy transfer. 



In this section, the criteria for the isoscalar nonlinear parameter 
values of the relativistic mean field model that cause instabilities at relatively high densities 
are investigated. 

We start from the energy density e of standard RMF models in symmetric nuclear matter 
(SNM) as, 

111 

£ = ^linear + ■^b2a^ + -630"^ - ^CsFq^, (1) 
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where eunear is the energy density of the hnear Walecka modeL Note that from now on, we 
denote parameter sets with nonhnear vector parameter C3=0 as S-RMF, and those with this 
parameter not equal to zero as V-RMF. The transverse (e^) and longitudinal (cl) dielectric 
functions at go = for SNM are [s^ 

er = l + 2rf^nT, 

eL = 1 + 2dsUs - 2d^Uv, (2) 
with the effective vector polarization defined as 

Uv = Hoo - 2dsUl, + 2dsUsUoo. (3) 
The explicit forms of the polarizations are 
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for the scalar polarization. 
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for the mixed scalar-vector polarization, and 
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for the longitudinal polarization, with Fermi momentum kp, nucleon effective mass M* = 
M - g^a, Fermi energy Ep= (kj, + M* '^fl'^, and E=(g74 + M* 2)V2, 
The longitudinal scalar meson propagator is given by 

qz _|_ ^* z 

while the vector meson longitudinal and transverse propagators are 

4 = 4= (9) 

gZ _|_ J^* Z 

where the a and a; mesons effective masses are given by 

^l'^ = =ml + 2b2a + 3b3a'^, 

= --^^ = ml + 3o,Vl (10) 

If we define m* ^ and ^ ^ as m* ^ and m* ^ at eT=0 and eL=0 simultaneously, then 
from Eqs. (jH [H] and [HD we have 



* 2 c 



(11) 



{Ut + Hoo) 

Because the presence of poles in the meson propagators is equivalent with er=0 and €l=0, 
it implies that at any nuclear matter density and momentum q, the stable regions are 
determined from the following criteria 



m* 2 ' 



m*, ^ 



< 1- (12) 



We observe that the m* ^ and m* ^ of every RMF parameter set strongly depend on 
the values of g„ and but are insensitive to the values of the parameters 62, and C3. 
In the contrary, m* ^ and m* ^ do not depend explicitly on g^r and g^^ values but they are 
very sensitive to the values of 62, ^3, and C3. Fortunately, we also observe that most of 
the standard RMF parameter sets in the literature's with acceptable predictions for nuclear 
matter and finite nuclei have similar g^r and g^j values. Let us take the average values of 
both coupling constants as the "standard" values of both parameters. We found that, for 
any parameter set, if the values of these (^o- and g^^ parameters are larger than the standard 



TABLE I: Isoscalar parameters for some test parameter sets. 



Parameter 


testl 


test2 tests test4 


test5 


9cr 


10.08 


10.18 10.02 10.02 


10.03 




12.62 


12.64 12.62 12.63 


12.62 


b2 


-10.53 


-7.10 -7.00 -6.80 


-6.83 




-28.00 


-3.80 0.00 3.80 


20.00 


C3 


0.29 


56.77 70.77 84.71 


159.70 



ones, its m* ^ ^ and m* ^ are higher and on the other hand, if these couphng constants are 
smaller, its m* ^ ^ and m*, ^ are lower. 

Based on the fact that almost all of the standard RMF parameter sets in the literature 
have similar values of and g^^ , we will use Eq. ( 1121) to determine a lower limit for the 62 , &3 , 
and C3 parameters of standard RMF models that are stable with respect to the particle-hole 
excitation modes at high densities. For the transverse mode, it is relatively straightforward to 
extract one parameter C3 from a critical effective omega meson mass, but for the longitudinal 
one, to pick up two parameters 62 and 63 from a critical effective sigma meson mass can 
generate many combinations because the correlation between 62 and 63 for each parameter 
set is different. It depends on many factors, like the values of other parameters, fitting 
procedure and strategy including the choice of observables for parameterization etc. To 
avoid this problem, we generate some test parameter sets with almost the same value of 
g^j and g^ with variation of 63 ( the 63 term dominates over the 62 term at high densities) 
while the 62 and C3 parameters are adjusted such that the standard SNM properties at 
saturation are fulfilled. The parameters of some exemplary test parameter sets can be seen 
in Table [H The SNM binding energy and pressure predicted by each test parameter in this 
table are displayed in Fig. [H while their ratioes of effective mass to nucleon mass and the 
compressibilities are shown in Fig. [2l It is clear from both pictures that they have similar 
nuclear matter properties at saturation (insets) but are different at high densities. The EOS 
of parameter sets with 63 > tends to be closer to the one of the experimental data from 

n ~ n 

Ref . [12| and the microscopic calculation of Ref . . 

In Fig. [3], the ratio of the critical effective meson masses to the effective meson masses is 
shown. The connection of the 63 parameter set with the appearance of instabilities is clearly 
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FIG. 1: The SNM binding energy (left panel) and pressure (right panel) calculated using test 
parameter sets in Table HI The Shaded region in the right panel corresponds to experimental data 
of Ref. [3], and the triangles to data calculated from the microscopic model of Ref. {l3]. The inset 
figure shows the low-density region at better resolution. 




FIG. 2: The SNM nucleon effective mass (left panel) and compressibility (right panel) calculated 
using the test parameter sets of Table HI The inset figure shows the low-density region at better 
resolution. 
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Pb/Po Pb/Po 



FIG. 3: m-CT ^ ^iid m* ^ '^/'"^^^ as functions of Pb/pq for various values of g and 63 as 

indicated. 
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FIG. 4: The onset of transverse instability (in the left panel) and effective omega meson mass (in 
the right panel) for some test parameter sets. The dot markers show their critical densities. 

captured. The explicit relation of the onset of transverse instability with effective omega 
meson mass can also be seen in Fig. H] and the one of the longitudinal mode is shown in 
Fig. [5l Fig. m shows the correlation between transverse instability and the effective mass of 
the omega meson which depends strongly on the vector nonlinear term. By comparing Fig. [6] 
with the values of C3 and 63 in Table [U it is clear that every parameter set with C3 > 71 does 
not show a transverse instability region at high density. To attain this value for C3, it seems 
that we need 63 close to or larger than zero. 

On the other hand, we can see from Fig. [5] that the longitudinal instability can be related 
to the effective sigma meson mass which strongly depends on scalar nonlinear terms. The 
longitudinal instability cannot fully disappear, but can be pushed to unphysical densities 
by taking 63 relatively large and positive, for example, 63 > 20, which it corresponds to 
the value of C3 > 160. In this case, the onset of instability can be pushed to pc > 23 po • 
For high-density applications (neutron stars, supernova matter, etc.) the maximum density 
is lOpo, with Po the SNM saturation density. This shows that for the purpose of these 
applications the parameter set with 63 > may be sufficient. For negative but relatively 
large absolute values of 63, an additional onset of longitudinal instability is generated. This 
appears due to the effective sigma meson mass becoming imaginary after reaching a certain 
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FIG. 5: The onset of longitudinal instability (in the left panel) and effective sigma meson mass (in 
the right panel) for some test parameter sets. The dot markers show their critical densities. 

critical density. This instability appears if 63 < —28. These results can be also directly 
seen from Fig. [31 Of course for arbitrary standard RMF parameter sets, its value can be 
relatively lower or higher, depending on how large the value of 62 is- Once again, we need to 
emphasize that the critical nonlinear parameters are obtained using "standard" and g^^ 
Therefore for parameter sets with and gc deviating appreciably from these, their critical 
nonlinear parameters can be smaller or larger. 

The last finding here is that the parameter sets which are stable at high densitiest are 
consistent with the experimental data of Ref. and the microscopic calculation of Ref. 
with respect to their SNM EOS. 



III. ANALYZING THE STABILITY OF SOME RELATIVISTIC MEAN FIELD 
MODELS 

Here we will use the criteria defined in the previous section to analyze more systematically 
the high density longitudinal and transverse instabilities of some RMF models. The onset of 



instabilities of some parameter sets have been reported in Ref. 



a 



but the discussion there 



was not quite robust. Therefore, in this work, we need to elucidate one point that was not 
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FIG. 6: The ratio between critical and saturation densities as a function of the 63 parameter for 
both modes. 

treated properly in Ref. 8|, namely that both instabilities are determined by the interplay 
between critical and actual effective meson masses (Eq. (fT2|) ). To comprehend our previous 
results js], now we also enlarge the number of parameter sets used. 



A. Standard relativistic mean field models 



First, we start with S-RMF. We select four parameter sets of this model as representa- 



tive, i.e., NLZ 



1^, NL3 



15|, NLSH 16| and NL2 [l^]. The isoscalar parameters of these 



parameter sets can be seen in Table [TTl It can be seen that and g,^ are of similar mag- 
nitude for NLZ, NL3, and NLSH, but the ones of NL2 are significantly smaller. Therefore, 
the parameter sets NLZ, NL3, and NLSH have similar a and uj critical effective masses 
while those of NL2 are smaller. Because NL2 has a and uj critical effective masses which 
are smaller than for the others and a positive value of the 63 parameter, the longitudinal 
instability of this parameter set disappears in the range of densities and q values used. On 
the other hand NLSH, having a quite large 62 parameter and a moderate 63 parameter, its 
sigma meson effective mass is larger compared to the ones of NLZ and NL3. Therefore, its 
instability regions are narrower and have no fluctuations. All parameter sets of S-RMF have 
a transverse instability. At relatively small g (< 4 fm~^), they start to appear in quite low 
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TABLE II: Isoscalar parts of some S-RMF model parameter sets. 
Parameter NLZ NL3 NLSH NL2 



10.05 10.22 10.44 9.11 

12.91 12.86 12.95 11.49 

62 -13.51 -36.12 -6.91 -2.30 

63 -40.22 -10.43 -15.83 13.78 
C3 



TABLE III: Isoscalar parts of some V-RMF model parameter sets. 



Parameter 


NLSVl 


PKl 


TMl 


NLSV2 


Z271 


FSUGold 


9cr 


10.12 


10.32 


10.03 


10.32 


7.03 


10.59 


9lu 


12.73 


13.01 


12.61 


12.88 


8.41 


14.30 


b2 


-9.24 


-8.17 


-7.23 


-6.86 


-5.43 


-4.28 


bs 


-15.39 


-10.00 


0.62 


0.37 


63.69 


49.93 


C3 


41.01 


55.65 


71.33 


72.39 


49.94 


418.39 



densities. These facts are clearly depicted in Fig. [71 

Next, we will discuss the V-RMF model. Here, we investigate the parameter sets NLSVl 
and NLSV2 of Ref. Q, TMl Q, PKl l2o|, Z271 [21], and FSUGold [4]. The isoscalar 
parameters of these parameter sets can be seen in Table lllli Due to the small values of 
the Qcr and g^^ coupling constants for the Z271 parameter set, it has relatively small a and 
u critical effective masses, so that even though its C3 < 71, the parameter set Z271 does 
not have a transverse instability region. On the other hand, for the TMl, NLSV2, and 
FSUGold parameter sets, their tranversal instabilities disappear due to the fact that these 
parameter sets have C3 > 71. In the V-RMF model, it is very clear reflected that, if 63 
becomes smaller then its longitudinal instability region becomes larger. This fact can be 
understood by comparing Fig. [8] with the 63 values in Table IIIII 

In general, the problem of the transverse part of S-RMF maybe can remedied by using 
very large bare omega meson masses but this still leaves the problem in the longitudinal part 
that should be fixed. Many parameter sets of the S-RMF model have negative but relatively 
large absolute value of 63. This fact leads to the situation that after reaching a certain 
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FIG. 7: m* ^ ^ /m*^ and m* ^ ^ /fntj' some S-RMF parametrizations in Table HIl as functions of 
PbIpqi with (7 varying as indicated. 
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FIG. 8: m* ^ '^/'m*^ and ^ ^/''T^t^ for some V-RMF parametrizations in Table Hill as functions 
of Pb/ Po and with q varying as indicated. 
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density their effective omega meson mass becomes imaginary. Tliis effect, as mentioned 
previously, appears as wild fluctuations in the critical effective sigma meson mass Q]. The 
presence of the C3 parameter in the V-RMF model leads to the point that some parameter 
sets of this model are free of transverse instability while for some the critical density for 
longitudinal instability can be pushed to an unphysical region. Thus the presence of C3 
larger than a certain value can indeed remedy both instabilities. 



B. Relativistic mean field model with mixing nonlinear terms 

In this subsection, we discuss in more detail than in our previous work Sj, the effect 
of the mixing nonlinear terms of some parameter sets of RMF model (E-RMF) that exist 
in the literature. The presence of these terms [sl] provides additional terms in the energy 
density e (Eq. ([T])) with —d2(TVQ — ^d^a'^VQ, and in the longitudinal dielectric function 
(Eq. ([2])) should be added with M^yUsv, where Usv = + '^d^v'^li + ^d^v'^s'^oo- The 
meson propagators therefore take the following forms 



9! 



^ g2 _^ ^* 2 _|_ ^2^(-g,2 _^ ^* 2^ 1' 



and now the scalar- vector coupling propagator is nonzero, i.e.. 



(13) 



SV (^2+^*2)(^2 + ^* 2)+^2^' ^ ) 



where the effective meson masses become 

A,^ = -{dh/dadVo) = 2d2Vo + 2d3aVo. (15) 
These terms also modify m* 2 in Eq. (ITT!) into 



m 



* 2 c 



n|, + nsnoo)(i + Ai) 



(nT + noo)(i + A2) 



(16) 



15 



while the corrections due to the mixing nonhnear terms are represented by 



= (n.n, - nU n.n„„) (^--"inL-n.n™]) 



^ 1 fx-n, 4:9ujgMTSo 



Ut + Hoo) V "'^ ^ A, 



cruj 



(17) 

with 

It is clear that if Ai = A2 = 0, these terms do not affect the critical sigma meson mass. If 
both A's are not equal to zero , the critical sigma meson mass reduction or increase depends 
on the sign and values of both parameters. On the other hand, the actual effective masses 
of sigma and omega mesons still depend explicitly on the nonlinear mixing terms. As the 
representatives of the E-RMF model, the parameter sets Gl and G2 from Ref. 22], as well as 



FSUZ003 and FSUZ006 from Ref. [3\ are used. The isoscalar parameters of these parameter 
sets can be seen in Table HVl The Ai, A2, m* ^ ^/m^ and Tnl^'^/rn^ of the representative 
parameter sets for the case q = 5fm~^ are shown in Fig. [9l In this case, the maximal 
values are Ai ~ 2.5% and A2 ~ 15% of the total contributions. Both are obtained using 
the Gl parameter set. This implies that the mixing nonlinear terms do not significantly 
affect the critical effective mass of the sigma meson. This can be seen from the lower left 
panel of Fig. [91 m* ^ '^/m^ of the Gl parameter set is significantly lower than those of the 
other parameter sets, merely because it produces significantly smaller g^j and g^^ compared 
to them (see Table II VI) . However, the significant effect of mixing nonlinear terms in the 
actual effective mass of the sigma meson can be seen in the right panel of Fig. [9l For the 
transverse mode, it is clear that G2 has an instability region because its C3 is critical, i.e., 
produces the maximal effective mass of the omega meson, equal to the critical value. If the 
parameter C3 is larger than this value like for the other E-RMF parameter sets, the instability 
disappears. For the longitudinal mode, the appearance of instability is merely determined 
by the effective sigma meson mass. The combination of 62, &3, and of each parameter set 
leads to mf^^ < m^^pg^^oos < ^*aG2 < "^^^Fsuzooe (right lower panel of Fig. [9]). Thus the 
critical density p^^ < p^suzoos ^ ^Gi ^ pFsuzoo6_ rjj^g large and positive value of 63 is still 
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FIG. 9: Ai (in left upper panel), A2 (in the right upper panel), m* ^ '^/m^ (in the left lower panel) 
and /f^t (i^ the right lower panel) for some E-RMF parameterizations and q = 5fm^^. 

the way to push the longitudinal instability to unphysical regions, but in the E-RMF model, 
the value of the critical 63 need not be too large, because its contribution can be partly 
substituted not only by that of 62 but also by that of ^3. These facts can be understood 
by comparing Fig. [TD] with nonlinear parameters in Table IIVI Another interesting finding is 
that if we observe the parameters of G2 in Table IIVI even its 63 is negative and its absolute 
value is relatively large, but other contributions from 62 and produce sufficient balancing 
contributions, avoiding this parameter set to have an imaginary effective sigma meson mass. 
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FIG. 10: m* ^ /"m*^ and m* ^ /nt^ of some E-RMF parameterizations of Table ITVl as functions 
of Pb/ Po and for some q values as indicated. 
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TABLE IV: Isoscalar parts of some E-RMF model parameter sets. 



Parameter 


Gl 


G2 


FSUZ003 


FSUZ006 


Qrr 


9.87 


10.49 


10.76 


11.02 


9ijj 


12.13 


12.76 


14.11 


14.66 


b2 


-15.09 


-24.89 


-9.72 


-4.65 


b3 


-47.69 


3.56 


21.80 


60.28 




86.41 


71.71 


198.25 


462.64 




-1.15 


-11.26 


-2.21 


-0.17 


d3 


-32.56 


4.19 


11.69 


1.38 



IV. CONCLUSION 

Based on the fact that the critical effective masses of the sigma and omega mesons (m* ^ 
and m* ^ ^) depend on the g^^ and couphng constants but are insensitive to the values 
of the 62, ^3, and C3 parameters, while the actual effective masses of the sigma and omega 
mesons (m* ^ and m* ^) have opposite properties, we derive criteria to estimate the values 
of isoscalar nonlinear terms of the standard RMF model that produce a stable EOS at high 
densities. The minimal requirement on the standard RMF model free from transverse mode 
instability at high densities is that the quartic vector nonlinear parameter should exist and 
have a value of C3 > 71. This value corresponds to the quartic scalar nonlinear parameter 
63 close to or larger than zero. The longitudinal mode instability cannot fully disappear but 
can be pushed to unphysical densities by taking 63 relatively large and positive, for example, 
^3 ^ 20, where in our chosen test parameter set, it corresponds to C3 > 160. In this case, the 
critical density for the instability region can be pushed to pc > 23 po • In usual applications, 
for examples neutron stars, supernova matter, etc., the density remains below lOpo? so that 
for these applications it seems sufficient to choose a parameter set with 63 slightly lower 
than 20 is sufficient. We have also found that the parameter set of the standard RMF model 
with a stable nuclear matter EOS at high density is consistent with the experimental data 
of Ref . and the microscopic calculation of Ref . [3] • 

The criteria has been used to systematically analyze the high density longitudinal and 
transverse modes instabilities of standard RMF models. We have shown that both parameter 
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sets unstable at high densities and stable ones exist. The reason behind this fact has been 
explained quantitatively and the crucial role of the isoscalar vector nonlinear term has also 
been clearly demonstrated. The effect of the mixing nonlinear terms of some parameter 
sets of the RMF model (E-RMF) has been also studied. Due to the contribution of these 
parameters to the actual effective masses of the sigma and omega mesons (m* ^ and m* ^) 
are changed, while the critical effective masses of sigma meson ( m* ^ are affected only 
weakly. The stabihty of both modes can be achieved more easily compared to V-RMF 
because the thresholds of the critical values of C3 and parameters of this model can be 
lower than those of V-RMF. 
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